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Abstract. We give an explicit formula for the effective partition function of a 
harmonically bound particle minimally coupled to a photon field in the dipole 
approximation. The effective partition function is shown to be the Laplace trans- 
form of a positive Borel measure, the effective measure of states. The absolutely 
continuous part of the latter allows for an analytic continuation, the singularities 
of which give rise to resonances. We give the precise location of these singularities, 
and show that they are well approximated by of first order poles with residues 
equal the multiplicities of the corresponding eigenspaces of the uncoupled quan- 
tum oscillator. Thus we obtain a complete analytic description of the natural line 
spectrum of the charged oscillator. 
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1. Introduction 

The standard model for a non-relativistic quantum particle interacting with the radia- 
tion field is the Pauli-Fierz model using minimal coupling. The Hamiltonian is 

H = ^-(p-eAf + V + H u (1.1) 

where Hf is the Hamiltonian of the free field, A is the transverse quantized field and e 
the electron charge, p is the particle momentum and V the particle potential. We do not 
introduce form factors. Over the last years, there has been much activity and significant 
progress in understanding the Pauli-Fierz model. Notable developments include the proof 
of existence for a ground state when the infrared cutoff is removed [H], a detailed spectral 
analysis of the Hamiltonian and the proof of the existence of resonaces [1] , and a proof of 
enhanced binding through the photon field [TO] . We refrain from giving a review of the by 
now vast literature on the subject, and instead refer to |18j and the bibliography therein. 

As is apparent already from the above selection of results, the majority of works on 
the subject investigates the spectral structure of the full system, particle and field. In 
our work, instead we regard the photon field as an environment in which the particle is 
embedded. A prominent example of this point of view is Feynman's polaron [8] by which 
our present investigations are inspired. The aim is to treat the coupled particle as a closed 
system, and to derive, for the particle alone, effective equations that capture the influence 
of the photon field. The particular effective quantity that we investigate (and, to our 
knowledge, indeed introduce) is the effective measure of states. 



We consider a version of (1.1) that is simplified in two ways: Firstly, we use the dipole 
approximation, which amounts to replacing A with A(0), and secondly we omit the self- 
interaction A 2 of the field. (The omission of the A 2 term is necessary only for a part of 
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our calculations, and we discuss the possibility of keeping it in the final remark of Section 
[3j) The particle potential V is taken to be harmonic, so that H is completely quadratic, 
and we add the well-known terms of mass and energy renormalization. In this setup we 
compute the ratio Z of the partition functions of the full system and of the free field. We 
do this by first suitably discretizing both systems, and then taking the continuous-mode 
limit and the ultraviolet limit. This procedure yields a finite result. We find the explicit 
expression 



P p —2p sin tp 
L27T 



T(ip e 



(1.2) 



with p ~ /3 inversely proportional to the temperature T and simp ~ e 2 determining the 



strength of the coupling. See Theorem |3.3| for details. 

From the ratio Z of the partition functions one gets an explicit expression for the 
difference of the respective free energies. Subtracting from this the free energy Fq of the 
particle leads to the excess free energy F ex = —(h) /3) ln(Z/Zo). The latter is particularly 
interesting in that it is supposed to be experimentally accessible. Referring to fjj we only 



mention that (1.2) confirms the well-known quadratic low-temperature behavior 



-7rak 2 T 2 (3mc 2 ) 



(1.3) 



of the excess free energy. 

In the present case, however, the significance of Z goes beyond that. As we will argue 



in Section 2.3 for general coupled systems, (1.2) is the effective partition function (or 



partition function for the particle subsystem). This means that Z replaces Zq when the 
small system is regarded as autonomous while retaining an effective influence of the field. 
The central fact is that, as we prove, Z is the Laplace transform of a positive Borel measure 
p. We interpret p as the effective measure of states for the charged oscillator. It has 
the same significance for the effective system that the measure of states po (whose Laplace 
transform is Zq, see (2.1) and (2.2)), has for the uncoupled oscillator. It determines the 
probabilities of energy measurements at the system in thermal equilibrium. 

It is worth remembering that in general the ratio of two partition functions is not the 
Laplace transform of some positive Borel measure. Most probably this is the case here 



as long as in the derivation of (1.2) the ultraviolet cutoff C is kept finite (even if the 
continuous- mode limit is already done). But a finite cutoff in not desirable in any case, 
as any result we obtain for the effective system would depend on the arbitrary parameter 
C. Thus it is crucial that the ultraviolet limit can be carried out, leading to an expression 
that is the Laplace transform of a positive measure. 

Formula (1.2) allows a detailed study of the effective measure of states. First of all, it 
turns out that p has one atom at the ground state frequency uj v , is supported on the half 
line [u^,oo[, and is otherwise absolutely continuous with respect to Lebesgue measure. 
The corresponding density 4> will be called effective density of states. We prove that <f> 
has an analytic continuation with singularities and cuts along the boundary of a cone with 
apex uj^ and opening angle 2cp. The conjugate pairs of singularities can be approximated 
by first order singularities in the interior of the cone, up to logarithmic corrections. The 
position and residue of these first order poles can be determined analytically. On the real 
line, they give rise to Lorentz profiles that replace the Dirac peaks present in p$. The 
mass of each Lorentz profile is equal to the mass of the corresponding Dirac peak, and its 
position agrees with predictions of QED time-dependent perturbation theory up to first 
order in the fine structure constant a. See Theorems |3 . 5 1 and |3 . 6| and the remarks following 
them for details. 
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From the above properties, a connection of the effective density of states with the theory 
of resonances appears obvious. The latter occur when a small quantum system is coupled 
to a quantum field or reservoir, at which point eigenvalues of the small system dissolve 
into the continuous spectrum, and stationary states change into metastable states. In the 
Pauli-Fierz model, resonances have been investigated extensively by Bach, Frohlich and 
Sigal in [lj. These authors employ the well-established complex dilation method. They 
study the spectrum of the full system, particle and field; resonances are then defined as 
the singularities of the analytic continuation of the resolvent to the second Riemann sheet. 
In contrast, we work in the context of open quantum systems, and therefore the effective 
small system does not even have a Hamiltonian. While there has been work on resonances 
and decoherence in open quantum systems, such as [2] and [TU [15], all authors seem to 
rely on the spectral structure of the full Hamiltonian. 

We define a resonance of an effective system in Definition [3j in terms of the complex 
structure of the effective density of states. For the harmonically bound charged particle, 
we show that this definition, and the concept of the effective density of states in general, 
produce physically reasonable results. Apart from facilitating explicit formulae, we do 
not expect the harmonic particle potential to be a vital part of the above picture; we 
thus conjecture that also for other confined charged systems, the effective density of states 
exists, and that its complex structure yields an appropriate description of the spectral 
lines by Lorentz profiles. A proof of this conjecture may well turn out to be challenging, 
and will probably require some insight into the connections of our results with the theory 
of P. 

Our paper is organized as follows. In Section[3j we present our model in detail and state 
the main rigorous results of this work, along with a short discussion. Proofs are given in 
Sections [4] and [5} The following Section [2] interprets our results in the framework of 
quantum statistical mechanics, and discusses their connection to the theory of resonances. 

Acknowledgments: We would like to thank Herbert Spohn for many useful discussions, and Marco 
Merkli for helpful comments. V.B. is supported by the EPSRC fellowship EP/D07181X/1. 

2. Partition function of the embedded system 

We regard the particle as embedded into the environment of the radiation field. An 
environment is characterized by the fact that it stays in thermal equilibrium if the joint 
system is in thermal equilibrium. Roughly speaking, the interaction has effects on the 
embedded system but does not change the environment. More precisely, when acting on 
the embedded system the environment carries out transitions between same states. 

A well-known example of a quantum system whose embedding in the environment de- 
termines decisively its properties is the polaron, i.e. an electron in an ionic crystal. The 
computations in [SJ Chapter 8] confirm that the electron moving with its accompanying 
distortion of the lattice behaves as a free particle but with an effective mass higher than 
that of the electron. In the same spirit we are going to treat the problem of a charged 
oscillator surrounded by its own radiation field. In this section we put h = k = 1. 

2.1. Energy distribution in thermal equilibrium. Let the positive operator H be 
the Hamiltonian of some quantum system. For the moment we assume that e~@ H is trace 
class. Classically, the partition function of the system is given by 



Z{p) = Tr e~ pH 



(2.1) 
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The inverse Laplace transform /i of Z, which we call the measure of states of the 
quantum system, is the weighted sum of Dirac measures 

M = ^iV i 5 A . ; (2.2) 
3 

where Xj are the eigenvalues of H and Nj are the corresponding finite multiplicities. Its 
physical relevance is due to the fact that the statistical uncertainty of states of the system 
in thermal equilibrium at the temperature T = 1//3 is described by the probability measure 
(x^ := j?py e~P('> fx. If E denotes the spectral measure of H, and A a Borel subset of K, 
then 

^ )(A) = zM Tr(e ^ HjB(A)) (2 ' 3) 

is the probability for a measurement of the energy to yield a value in A. It is this 
interpretation that we will retain valid also for the effective entities. 

2.2. Spectral discretization. From the basic principle of statistical mechanics, once 
the partition function is known, all thermodynamic properties can be found. In many 



interesting cases, however, including our model Hamiltonian (3.1), e~@ H is not trace 
class. This raises the problem of how to attribute a partition function to the oscillator 
embedded in the photon field. The way we solve this problem uses a slight generalization 
of methods from the theory of random Schrodinger operators (cf. [12]). The main idea is 
to approximate H by operators H n with discrete spectrum. We define 

Definition 1. Let H be a self-adjoint operator in a Hilbert space H, and let (P n ) be a 
family of orthogonal projections on H. Put H n = P n HP n . We say that H n is a spectral 
discretization of PL ( associated to P n ) if 

(i) P n D(H) C D(H) 

(ii) lim ri _ 5 . OC) P n = 1 in the strong topology. 

(iii) The spectrum of H n acting in P n H consists of eigenvalues with finite multiplicity 
such that e~P Hn is trace class for all n. 

2.3. Separation of a subsystem. Let us now consider a general model consisting of two 
interacting quantum subsystems with respective state spaces Hi and Hi. For convenience, 
the system corresponding to Hi will be called the small system, the other one the large 
system; but note that there is actually no assumption on the dimension of Hi and Hi. 
The Hamiltonian of the joint system, acting on Hi ® H2 is 

H = Hi ® 1 + 1 <g> H 2 + Hi, (2.4) 

where Hi describes the interaction. Again we assume for H the existence of Tr e~@ H , 
and also for Hi and Hi. Let the large system be in thermal equilibrium. Then the 
influence of it onto the small system can be described in a statistical way by averaging 
the states of the large system. Accordingly, the density operator attributed to the small 
system is W(/3) := (Tr e~@ H ' 2 ) _1 Tr 2 e~^ H , where Tr2 denotes the partial trace with 
respect to the second factor. In the case H\ = of non-interacting systems, W(J3) equals 
Q—pHi _ Moreover, if the joint system is in the state (Tr e~@ H )~ l e~@ H , any prediction 
on measurements which concern only the small system is given by the thus uniquely 
determined state (Tr W{P))~ X W{P) = (Tr e _/3H ) _1 Tr 2 e _/3H of the small system. This 
is due to the fact that Tr (Tr i^~^ H E\) = Tr(e~^ Ei ® 1) holds for every orthogonal 
projection E\ on Hi and uniquely determines Tr2- 
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Consequently, we regard Z{j3) := TrM^(/3) as the partition function attributed to the 
small system. It satisfies 

Tr e'P 11 

m = (2-5) 

and we will interpret its inverse Laplace transform as the effective measure of states (cf. 



2.2) for the small system when in contact with the large system. Of course, for this 
interpretation to make perfectly sense, Z would need to be the Laplace transform of a 
positive Borel measure. This is not true in general, but we will find that it does hold in 
our model when performing the following limiting process. 



2.4. Infinitely large environment. The large system in (1.1) is described by H2 = H{, 
which has absolutely conti nuo us spectrum on [0,oo]. Thus e~^ H2 is not trace class, and 



neither is e ^ , whence (2.5) is not available. We consider spectral discretizations H n 



and i?2,n of H and H2, and define 

Z M = xT^r- (2 ' 6) 

In the case we study lim ri _ Sl00 Z n exists. After removing an ultraviolet cutoff we get the 
partition function Z, which is the Laplace transform of a positive measure \x with support 
in [0,oo[. It is called the effective measure of states of the embedded system. As in 
Section 2.1 we introduce the probability measures 

u (ffl - _J_ e -P(-) a (2 7) 

for f3 > 0, which are fundamental in that they determine the probabilities of energy 
measurements: if the charged oscillator is in thermal equilibrium at the temperature -g, 

then n^\A) is the probability for a measurement of the energy to yield a value in A. 

One can read off Z and \i the same amount of information about the spectrum of the 
embedded system as in the case of a trace class operator. We will see that the embedding 
into its own radiation field changes the behaviour of the oscillator qualitatively: instead of 
a purely discrete spectrum, we now obtain apart from the stable ground state an absolutely 
continuous spectrum on the positive half axis. 

2.5. Resonances. An interesting property of the effective measure of states [i of the 
present system is that the analytic continuation of its density function <j) has conjugate 
pairs of first order singularities, which for small coupling are very close to the real line, 



cf. Theorems 3.5, 3.6 They manifest themselves as Lorentz profiles (Breit-Wigner reso- 
nance shapes) in the effective density of states. There is a connection with the theory of 
resonances |11| . which we will elucidate here. 

Generally speaking, bound states, which are perturbed, give rise to resonances. As 
shown in [11] , [1] , in case of a bounded electron coupled to a photon field these are related 
to singularities p of the resolvent of the Hamiltonian H off the real axis coming from the 
eigenvalues of the uncoupled electron system Hq. In accordance with the definition of a 
resonance given in |16| XII. 6] this means that there are p 6 C with Im p < and a 
dense set T> of state vectors u for which the matrix elements R> u \z) = (u, (H — z) -1 -u) 
and R < ^ 1 \z) = (u,(Hq — z)" 1 ^) have an analytic continuation from the upper complex 
half-plane across the positive real axis into the lower complex half-plane such that the 
points p are singular points of the former but regular for the latter. 
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There is an equivalent formulation in terms of the scalar spectral measures. Let E be 
the spectral measure of H. For given state vector u, denote by the scalar spectral 
measure on [0, oo[ given by fi^ u \A) := (u, E(A)u) . Let V be the maximal open set 
on which the distribution function w of is real analytic. We consider the complex 
analytic continuations f"' of w'\V on symmetric domains with respect to the real axis so 
that 0(«)(z) = <f>W(z) holds. Let cp^ denote the respective object for Hq. 

Definition 2. Let T> be a dense set of state vectors andp = ii-i^6C with u) > 0, 7 > 0. 
Then p is a resonance if for every u € T> there exist complex analytic continuations cf)( u ' 

(u) 

and (J)q such that p is a singular point of the former and a regular point of the latter. 

In the Appendix it is shown that the two definitions of a resonance are equivalent. More 
precisely, R^ is analytic on C \ [0, oof and its analytic continuation across the positive 
real axis to the second Riemann sheet equals + 27ri0 ( - u \ The same holds true for the 
respective objects for Hq. 

The singularities of occur in complex conjugate pairs. They manifest themselves 
as Lorentz profiles in (jy^) along the real axis. For an illustration assume the simplest case 
that the resonance p is a pole of first order and that there are no other singularities than 
p and its conjugate p in some open rectangle U, which is symmetric with respect to the 
real axis and the axis through p, p. Then there is a holomorphic function <p( u > on U such 

that 4>( u \z) = 7^pj7^~p) f° r z £ U \ {p,p}. Hence along the real axis not far from Cj one 
has 

L(U)(. ^ 7/2 



4> w {u) ~ constant y- .., . . 

which is a Lorentz profile. A quasi energy eigenstate u, which in the energy represen- 
tation shows a narrow frequency band A around u>, decays exponentially. Indeed, the 
matrix element (u, e~ ltH u) is given by the Fourier transform pS u '{t) = J e -1 ^ dfj,( u '(u)) of 
the scalar spectral measure, which in the present case becomes 

= f e~ ituJ <j){Lu)du ~ constant e"^* e~* 7/2 . 
J A 

The matrix element R^ (z) of the resolvent equals the Stieltjes transform of the scalar 
spectral measure 



uj — z 



The formula 



1 f°° 

jj,(t) = — lim e _;.o+ — / e -1 *" 1m p,(ui + ie)dw, 

7T J -co 



which converts the Stieltjes transform of any finite Borel measure fi on [0, oof into its 
Fourier transform, allows to study the decay of a state u related to a resonance starting 
from J2("). Cf. pj 1.32] and [TI]. 

To summarize so far, we have seen that it is quite natural (and equivalent) to determine 
the resonances by the analytic continuation of the densities of scalar spectral measures 
of the Hamiltonian rather than by the analytic continuation of matrix elements of its 
resolvent. Now the question is how to apply these considerations to an effective system 
for which there is no Hamiltonian. Note that once the effective measure of states [i is 
determined, the probability distributions fiW of the energy are at our disposal for every 
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temperature 1/(3 (see Section 2.4). These are the primary entities in case of an effective 
system. It appears natural that they have to replace the scalar spectral measures of the 
Hamiltonian. Hence common singular points (in the lower complex half-plane) of the 
analytic continuations <f>^' of the respective densities are attributed to resonances. It 



is immediate from (2.7) that the former are exactly the singular points of the analytic 
continuation (j) of the density regarding the effective measure of states [i. Thus we are led 
to the following 

Definition 3. Resonances of an effective system are the singularities (in the lower complex 
half-plane) of the analytic continuation of the effective density of states. 

We have seen that resonances according Definition [3] give rise to Lorentz profiles along 



the real axis representing the natural lines of the system. The results in Theorems 3.5 and 



3.6 then confirm that these resonances actually occur for the charged harmonic oscillator. 



What is more, they furnish a physically verifiable picture of the spectrum in its entirety. 



3. Model and main results 

3.1. Particle coupled to a photon field. Our starting point is the minimal coupling 
without a form factor (cf. [U II.D.l (D.l)]), in the dipole approximation of a charged parti- 
cle interacting with a photon field and in a quadratic external potential. The Hamiltonian 
is given by 

H = (H p + H T )®1 + 1®H { + Hi (3.1) 
acting in L 2 (M 3 ) ® J 7 ® 2 . Let us explain the symbols above. 

acting in L 2 (M 3 ), is the energy of the harmonically bound particle, m is the particle mass, 
and r\ is the frequency associated to the the harmonic external potential. Note that we 
represent the particle in momentum space, which will later turn out to be convenient. 

Hi = H fi i <8> 1 -(- 1 <g> H fj2 and H f)(J = hj w(k)a*(k, cr)a(k, a) dk, a = 1,2 

is the energy of the free field. Each Hf t<J , a = 1,2, acts in the symmetric Fock space 

? = 0n=o- F(n) > where • F(n) is the set of a11 / G L 2 (R 3n ) such that f(h,...,k n ) = 
/(&7r(i)) • • • j ^7r(n)) f° r an permutations it, and all kj G M 3 . We write an element F of the 
Fock space as F = J2n=o F(n) with E^°=o ll^ (n) l| 2 < °o for F^ E F^ n \ We will need the 
following well-known fact: if {f m : m G N} is a orthonormal basis of L 2 (M 3 ) or a subspace 
X thereof, then the system 

{/mi® • ■•®/m„ : "il, • • • ,m n G N, n G N} 

is a orthonormal basis of the Fock space over L 2 (M 3 ) or X, respectively. Above, <S> denotes 
the symmetric tensor product. 

By definition, Hf CT acts on J 7 ^™) as 

(h f uj(k)a*(k,a)a(k,a)dkF ( - n A (k l ,...,k n ) = Y,hu}{k j )F ( - n \k 1 ,...,k n ), 

v J ' j= i 
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where oj{k) = c\k\ is the photon dispersion relation with c the velocity of light. The 
interaction between the particle and the field is given by 



Hi = Hi tl ®l+l®Hi :2 and H^ a = ~ J Xc(k)J (q-u a (k)) (a* (k, a)+a(k, a)jdk, 

with an ultraviolet cutoff implemented through Xc = l{|fc|<c}- It is wen known (see e.g. 
[3]) that Hi is infinitesimally form bounded with respect to H p ® 1 + 1 ® H{ , so that H is 
essentially self-adjoint on D(H p ® 1 + 1 ® Hf). 

The electron charge e determines the strength of the coupling, and the transverse vector 
fields u a , such that u±(k), U2(k), k/\k\ are orthonormal, account for photon polarization. 
Hi acts on L 2 (dq) as a multiplication operator, and on F the creation and annihilation 
operators a*(g) = J g(k)a* (k) dk and a(g) = J g(k)a{k) dk act in the following way: The 
creation operator takes F^ into 

jr(n+i) through 

1 71+1 

(a*(g)F)^ +1 \k 1 , k n+1 ) = -=== V g{k i )F^\k x , . . . , fe, . . . , k n+ i), 

where \j means that the argument kj is omitted. The annihilation operator takes F^ 
into J 7 ^ -1 ) through 

(a(g)F)^~ 1 \k l , fc n _x) = g(k)F^ n \k, fci, . . . , fc n _i)dA;. 

Finally, 



2vr 2 m J u 2 (k) y ^ 4vr 2 



' (eVHCg 2 -^(eVHln(l + -) 



07rcr 7rc 7/ 

is a renormalization acting nontrivially only on L 2 (dg). The first term yields the stan- 
dard mass renormalization accounting for the increased mass of the particle due to the 
coupling to the photon field, and the second term is the well-known energy renormaliza- 
tion, which compensates the negative shift of the energy scale also due to the interaction 
with the photons. Including these two terms will keep finite results when we remove the 
ultraviolet cutoff by taking C — > oo. 

We now give spectral discretizations for H and H{. For N 6 N let Aj, j £ N, be the 
half open cubes with volume V = N~ 3 and vertices in the lattice (^Z) 3 C M. 3 . We define 
the one dimensional orthogonal projections 

P N (j) : L 2 (M 3 ) -+ L 2 (R 3 ), /^(hjj)^ with h s := ^U r (3.2) 



Let Jjv := {j G N : n B(0,N) ± 0}, where B(0,N) is the centered ball with radius 
N. Since P N {j) P N (f) = for j j' , 

Pn ■= E p ^ 

J'GJjv 

is a finite dimensional orthogonal projection. Hence so is .Pat, the second quantisation of 
Pn, acting on F^ as 

Pn /i® • • ■ <§>/„ = (Pjv/i)® • • • MPnU). (3.3) 
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We then define 

#f,N := (Pn ® Pjv) Hf (P N ® Pat), iljv := (1 ® -Pjv ® Ptv) H(l(g)P N ® P N ). 

Proposition 3.1. ^ and Hn are spectral discretizations of Hf and H. 

The proof will be given in Section |4j 

3.2. Main results. Our main results are about the partition function Z(JS) of the oscil- 
lator embedded in the photon field, and about its inverse Laplace transform, the effective 
measure of states. It is convenient to introduce dimensionless entities. The strength of 
the coupling is given by an angle ip such that 



srntp :- 



3ucf 



a 



3mc 3 



holds, where a 



he 



is the Sommerfeld fine structure constant and locf 



mc" 

h 



(3.4) 
the Comp- 

ton frequency. As to the cutoff C we introduce 7 := cC/n, and (5 = ^ is replaced by the 
variable p := 

Theorem 3.2. Tr e~f Hjv and Tr e~^ f > N exist for each N G N and each f3 > 0. More- 



over, Z{J5\ 7) 



lim7v^oo(Tr e""*^ )/(Tr e~h H t^ ) exists for each (3 > 0, and 



2irp e 



-Ip ln(l+7) sin ip 



n 



p z 4 p , , . ( p 

1 + -s- H — sin(</?J arctan 7- 

l z IT I \ I 



(3.5) 



Formula (3.5) already appears in [4], where it is derived non-rigorously using path inte- 
grals. We present a conceptually simpler, and mathematically rigorous, proof of Theorem 
531 in Sections ISTTl-IOl 



Removing the ultraviolet cutoff in (3.5) results in a significantly simpler expression. 
Theorem 3.3. For all ft > 0, Z(/3) := linXy-^oo Z(/3;7) exists, namely 



Z{P) 



P -2pln(p)siny Ipfj/, g^i^)! 2 
2vr 1 VP ;| 



(3.6) 



Moreover, (3 1— > Z(/3) is i/ie Laplace transform of a positive measure. 



The statements of Theorem 3.3 also appear in [3]. But the derivation of (3.6) is not fully 
rigorous there and the proof of the final statement contains a small error. Both results 



are improved in Sections |5.4| and 6.1, respectively. 

We now turn to the effective measure of states p v of the charged oscillator, i.e. the 
Laplace inverse of Z v . We display the dependence on the strength of the coupling by the 



suffix <p € [0, ?]. For <p = 0, we are in the case of zero coupling, Zq{(3) 
the partition function of the oscillator, and 



(2 sinh 



2 / 



-3 ; 



is 



Mo 



j=0 



j 2 2 



5 (j+3/2)»r 



(3.7) 



By (2.7), ^ 



ature T 



p v is the probability distribution of the energy at the temper- 



An easy first step is the following convergence result. 



Proposition 3.4. As (p 

the vague topology. 



0, convergence pip 



p^ and p v 



Pq holds in the sense of 
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Proof. One finds (i^\t) = Z^{fi + it)/Z^((3) for the Fourier transform of iffl at t € R. 

Since ^(/3 + it)/Z,p(f3) — >• Zq((3 + H)/Zo(/3) = Ao (*) as (/?—)■ 0, the vague convergence of 

((J^p^)ip follows from the continuity theorem. This implies the vague convergence of (fJ,<p)(p 
simply noting that if ip is a continuous function on R with compact support then so is 
exp(/3 (•))</>. □ 

The foregoing result implies that the probability mass of fi^ concentrates near the lines 
at (j + 3/2)77 and vanishes in between as the coupling strength tends to zero. This is 
exactly what is expected. In the following we study the structure of /j^ more closely. We 
define 

ujy := - (sirup + - ip) cos ip)r], (3.8) 

7T V 2 J 

and write 6 U for the Dirac measure at w^, and c^A 1 for the Lebesgue measure on R with 
density (j). 

Theorem 3.5. 

i) There exists a function <p : R — > R being zero for uj < lu v and positive continuous 
for uj > uj v with (j){oo^) = — sirup such that 

ii) <p is real analytic for oj > uj v , and extends to an analytic function on C up to 
singularities at 

pj =u tp -\- jr/ e~ llfi and pj = u) v + jr] e 1(p , j G Z \ {0}, 

and cuts along z = ui^ + s e ±1</3 for s£R, \s\ ^ n. 

So the zero-point frequency |r/ of the oscillator is shifted down to oj v when the oscillator 
couples to the radiation field, and (p is the effective density of states of the charged 



oscillator. Part i) above is already shown in [1], but in Section 6.2 an elementary proof is 



given. According to i) there is a stable ground state at the energy hui^, while at the same 
place the absolutely continuous part of the spectrum jumps to ^ sin^. Hence the ground 
state is no longer isolated. It also implies that there is no other stable and no singular part 
of the spectrum. Part ii) goes beyond these qualitative statements, and already shows the 
existence of resonances (cf. Definition^) at pj, j £ N. Its proof will be given in Section 



6.2 However, it still contains no information on the nature of the singularities; in order to 
observe unadulterated Lorentz profiles on the real line, these need to be simple poles. 
This, and more, is true. To formulate the corresponding result, let us define 

:= 2m(z- Pj ) + 2vri(z - pj) ' 
for z S C. When restricted to the real line, 

1 jr} sin ip 



ir (oj — bj^ — jr] cos tp) 2 + (jr] sin <p) 2 
is a Lorentz profile with total mass J £j du = 1 . 
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Theorem 3.6. 



Let N eN and z 

N 



+ se' x for < s < Nrj and \x\ < 



Then 



E 

3 



i + 2 

2 



£,•(» + /iat(z) i/|x| < <f, 



^(-iy( 3 \ 3 ( z ) + h N (z) if\x\><p. 



Above, hjy and hpj are analytic up to singularities atpj, pj and cuts along z = w^ + se ±¥ , 
s ^ rj. Moreover, there exist a constant C depending on N such that hjy (and similarly 
h N ) satisfies \h N {z)\ < C(l + | hx(<p - \ X \)\ vN ) for \ X \ < V- 



The proof of this theorem is given in Section 6.3 The theorem shows that inside the 
cone |x| < ¥?, the singularities of cj) look like resonance poles of first order, but their overall 
structure is more complicated as is apparent from the behaviour for \x\ > y>- Nevertheless, 
they lead to Lorentz profiles ( 5 ^ 2 )^j on the real axis. The total mass of the j'-th Lorentz 
profile thus corresponds precisely to the multiplicity of the j-th eigenvalue of the uncoupled 
harmonic oscillator, i.e. the mass of the j'-th delta peak of no, for all <p. Moreover, with 
Proposition 3.4 it follows that duo converges vaguely to zero as <p — > 0. 

To summarize, the following complete picture of fi v emerged: There is a stable ground 
state (Dirac peak of strength one) at hoJ v . Above, there is absolutely continuous spectrum, 
consisting of Lorentz profiles with widths 

2a ( Hri \ 

7j := 2jr)sm<p = jn— [ ) , (3.9) 



centered at 



a 
18 



Hr] 



+ 



(3.10) 



and with total mass equal to that of the corresponding unperturbed state. 



150 



50 - 





b) M/2 



Figure 1 . a) shows the effective density of states, computed by summing 



j+l/2 



up the first 15 terms of (6.4), for (p = 1/100. All terms of (6.4) are posi- 



tive, so the numerical approximation is very good, b) shows the difference 
between the effective density of states and its approximation by the appro- 
priate Lorentz-profiles (J + + 2)£j/2, in the intervals [j — 1/2, j + 1/2], 
for j = 1, 2, 3, 4. In both graphs, the units of energy on the x-axis are cho- 
sen such that the eigenvalues of the uncoupled harmonic oscillator occur 
at the elements of Nq. 
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Figure[T]illustrates the excellent approximation to the effective density of states obtained 
by the above description, even for the relatively large value tp = 1/100; the physical value 
for ip is about a~ 3 « 0, 5 x 10~ 6 . 

With respect to a frequency scale starting at uj v the peaks at ojj for j = 1,2,3,... 
represent the natural lines of the charged oscillator due to transitions to the ground 



state (Lyman series). According to (3.10) the distance of two subsequent levels is rjcostp, 
which is smaller than 77 in the unperturbed case. The radiative shift (Lamb shift) 77— r/ cos tp 
of the first excited state is of second order in a. 

For an observation of the natural lines the charged oscillator may be brought into contact 
with black body radiation at the temperature T = A. This causes additional negative 
shifts of the excited states depending on T which have to be taken into account. They 
are determined by the positions of the peaks of /4? . For IT > (ft/k)7j one readily finds 
-07,74 < A Wi < -^/8. 

Our results are in accordance with calculations within QED time-dependent perturba- 
tion theory up to first order in the fine-structure constant a (cf. e.g. [5j sec. 5]). The 
calculations on the radiative cascade of a harmonic oscillator in [7, Exercises 15] confirm 
also that the shifted levels of the oscillator remain equidistant. 

Remark. In our analysis, we have left out the self-interaction term of the field (A 2 )- 
term). At least Theorem 3.2 can be achieved without this simplification. Actually, in the 
Goppert-Mayer description this term is absorbed by the coupling, compare [6j IV (B.34)] 
or (TBI (13.124)] and [13 V.15 eq. (15.1), (15.2)]. One is left again with an oscillator 
model which can be treated in the same way as we do here. The analogous formula to 



(3.5), without adding any extra renormalization term, is 



1=1 



sin ip 7 



I ( p 

— arctan I 7- 
p V I 



Again the limit 7—7-00 does not exist, but this time we do not know how to suitably 
renormalize this system. 

On the other hand, from the well-tested perturbation theory in QED one knows that 
for many problems concerning the interaction of bound electrons with radiation (e.g. 
spontaneous and stimulated emission and absorption) the linear part — epA/m of the 
interaction Hamiltonian yields the main effects (see e.g. [13]). 



4. Discretisation of the photon field 



Let us now give the proof of Proposition 3.1 and at the same time prepare the proof of 



Theorem 3.2 In effect, what we give is the reverse procedure to the one that, in the early 
days of QED, led from the description of the photon field as a collection of independent 
harmonic oscillators to the Fock space description. So in principle, what we are going to 
present is well known. However, we are not aware of any place where it is done carefully 
and in a mathematically satisfactory way, and thus present it in some detail. 



Let us start with some properties of the projections Pjy introduced in (3.3). For locally 

p3 



integrable / 



C, and given j £ Jjv, we write 
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for the average value of / on Aj. Thus, if / G I? we have 

Pn/(x)= fM x )> x£R3 - 

Lemma 4.1. Assume f G L 2 (R 3 ), F G F . Then we have: 

(i) As N -> oo, P N f -> f in L 2 (R 3 ), and P N F -> F in F '. 

(ii) If f is continuous on an open set B, then P]\rf(x) — > f(x) for every x G B. 

Proof. As to (ii), assume without restriction that / is real- valued. Choose xo G B. For all 
N large enough, we have Aj C B for the unique jo such that xq G Aj . By continuity there 
are x N , x' N in A io satisfying f(x N ) = inf /| s and f(x' N ) = sup/| s . Then P N f(x ) = 

fjo e [f( x N),f(x' N )], whence the assertion. 

Now (i) is proved in the standard fashion. Choose / to be the indicator function of a 
bounded measurable set B at first. Then for each e > there is a continuous functions g 
with compact support such that \\g — 1b\\l2 < £■ Now 

\\Pn1b-1b\\ < \\Pn(9-Ib)\\ + \\Pn9-9\\ + \\9-1b\\. 

The middle term converges to zero by (ii) and dominated convergence, and so the left hand 
side is bounded by 2s for all large N. This proves the claim for indicators of bounded 
measurable sets. As the latter are total in L 2 , this proves also the general result about 
(Pn)n- Finally, the result about (Pn)n follows from the fact that second quantisation 
preserves strong convergence. □ 

It can be checked easily that PnJ 7 = J'iPNL 2 ^ 3 )), the latter being the Fock space over 
the image of Pn- Since {hj : j G Jn} is a ortonormal basis of PnL 2 (W 3 ), we conclude that 
{hfa® . . . ®hj n : jk G Jn, n G No} is a orthonormal basis of PnF- These basis elements 
are eigenfunctions of PnH{ j(T Pn, as 

P N H { >CT h h ® . . . ®h jn = (^ h Yl Q ii j h h ® ■ • ■ ®h jn . 

This allows us to transfer the operator Hf^ unitarily to an L 2 space. For n G N we write 
for any positive mass parameter fx 

ipn{x) = 

for the n-th Hermite function on R. For j G Jn and j := (ji, . . . ,j n ) G J N , n G N let 
n(j,j) denote the number of occurrences of the value j in j; for n = set j := 0, := 
1 G C and n(j, 0) := for all j G J N - We put 

hj := h h ® . . . ®h jn G F^ and ^ := (g) G L 2 (R Jn ). 

Due to symmetrisation, hj actually only depends on the occupation numbers n(j,j), j G 
Jn, and thus it is easy to check that the asignment 9(hj) := ipj defines an isomorphism 
from P N F onto L 2 (R Jn ), and that 

ep N H, a e-i = Y: (-fol 3 + \*)xl 3 - *f) , (4.i) 
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where x\j and x 2 j are independent variables. Due to the relation 2xH n (x) = H n+ \{x) + 
2rt-ff n _i(ic) we have 

Vn + iip n+1 (x) + V^V'n-i(^) = J-^-^n(x), 
and using the definitions of a* (g) and a(g), it is tedious but straightforward to check that 

0P N Hi, a 6~ l = J^2~^ E V^^i ' ^ Xu 3 witn y° := u <tXc/Vu- 



J'GJjv 

Thus we find that is unitarily equivalent to 



acting in L 2 (R 3+2Jn ). Here we discretized also the renormalisztion terms 



(4.2) 



^n=E E^fe^-,) 2 - E™; 



displaying the contribution to the effective oscillator mass and energy shift by each single 
photon. 

Since Hn is the Hamiltonian of a finite dimensional quantum oscillator, its spectrum is 



discrete and the eigenvalues are of finite multiplicity. Using Lemma 4.1, it is immediate 



to check (i) and (ii) of Definition [I] for 1 ® Pn &> -Pn> and Proposition 3.1 is proved. 

5. The partition function 



In this section we will prove Theorems 3.2 and 3.3 We start with a general, useful 
observation. 

5.1. A trace formula for a quantum oscillator. Consider an operator 

n _ %2 



2mi 1 2 

i=l i,i'=l 

with strictly positive real n x n-matrix A = (An'). H acts in L 2 (M n ), dt are derivative 
operators with respect to X{, and X{ are multiplication operators with Xi- We assume 
that the mass parameters rrii are all positive, and we define the diagonal mass matrix 
M := diag(mi, . . . ,m n ). 

Lemma 5.1. Let Ai,... 5 A n be the eigenvalues of M~ 1 / 2 AM~ l l 2 . Then 

Tr e-* H = f[ ^sinh^v^)) ' = det ^2 sinh |^(M _1 ' 2 AM 4 ' 2 )^j ' = 

r n det{M-y 2 AM-^ 2 )-^ f[det (l n + {^fM'^AM-^A * . 

i=i ^ 71 ' 
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Proof. We are going to show that H is unitarily equivalent to 



frf \2m 

i=i v 



~2~ * 



(5.2) 



with mo > any mass parameter. Then by Mehler's formula the first equality holds. The 
second equality is obvious. The last equality holds true due to the factorization 



sinh z 



n(»+< 



1=1 



valid for all z G C. To show (5.2), we note that any invertible real n x n-matrix S gives 
rise to a unitary transformation U on L 2 (M n ) by Uf := \detS\% f ° S. We choose S in the 
following way: Since M~ 1 / 2 AM~ 1 / 2 is real symmetric, there exists an orthogonal matrix O 
such that OM~\AM-\o T = diag(Ai, . . . , A n ). We put S := 0(^M)i Then S satisfies 



SM^S 7 * 



m 



-In, 



S-^AS' 1 



m diag(Ai, . . . , A n ), 



and from these relations H = U 1 H U follows straightforwardly. 



□ 



5.2. The Partition Functions for H{ jv and H^. We will now use Lemma 5.1 in order 
to compute the partition functions for the discretized systems. Note that from ( |4.T| ) and 
Lemma 5.1 we immediately get 

-2 



Tr e 



J12 



n ( 2sinh (f w i' 



(5.3) 




Now by (4.2), Hjy has the form (j 5 . 1 p up to the constant term 

(Xc)j \ 



u>i + V 



e 2 i] 2 



4ir 2 m Lo 2 (ujj + rj) 



Precisely, the mass matrix M is the diagonal matrix of size 2| J^ \ + 3, with as the first 
three diagonal elements, and /i as the remaining ones. A is a block matrix of the form 



(5.4) 



with B = -jri3 (/j denotes the j-dimensional unit matrix), Sj = [iuj 2 - for 1 j ^ \ Jn\i an d 





f B 


B 1 


B 2 


' ' ' B \j N \ \ 






s\h 










Bl 












\ B \Jn\ 







S I^I J 2/ 



= (vy,uai) G M 2x3 , with v. 



2tt 2 m ' 



foJjyaj for <t = 1,2. It is clear that both of 
the matrices whose determinant appears in the second line of Lemma 5.1 are of the above 
form, too. We thus need a formula for determinants of matrices of the type (5.4). 



Lemma 5.2. Let Q be a real symmetric {2L + 3) -dimensional matrix of the form (5.4), 
where B is 3- dimensional and all Sj ^ 0. Then 



detQ 



sj det B 
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Proof. We write Bl = (a,b), where a and b are 3- dimensional column vectors. We 
assume for the moment that a and b are linearly independent. At the end this assumption 
can be dropped by continuity. Define a 3 x 3-matrix X by X T := (x,y,z) with x := 
a x b, y := b x x and z := x x a. Then XBl = \a x b\ 2 (e2,es) with I3 = (ei,e2,es) and 
X -1 = \a x b\~ 2 (x, a,b) hold. The following operations do not change detQ. 

First pass to the equivalent matrix diag(X, I2L) Q diag(X , Then the last block 
in the first block line becomes zero adding —\a x b\ 2 /si times the last two lines to the 
second and third line of the matrix. Expanding now the determinant along the last two 
columns one obtains the factor s\ times the determinant of a (2L + l)-dimensional matrix 
Q' . One finds that diag(X _1 , 12L-2) Q' diag(X, I2L-2) arises from Q by canceling the last 
block line and last block column and replacing B with B — j^B^Bj^. Iterating these 
operations the result follows. □ 



Thus by elementary computations, Lemma 5.1 and Lemma 5.2 yield the following for- 



mula for the ratio of the partition functions of the discretised systems: 



Z N {fi) :-- 



Tr e" 



Tr e"R fff - N 



/9E 



2 2 
y e V 

J£ J N 4w 2 m ~E/Z~< 
3 1 



(Xc)j 



i+V) 



J\det 



e 2 r] 2 



1=1 



2ir 2 



V 



m 



E 



VajV, 



T 



(5.5) 



E ^ , . , 

j£j N a=l,2 Jy 3 I 



with 



Va ■= Xc 



Ho 



UJ 



W 



5.3. The continuum limit. We will now show that in formula (5.5), the limit N — > 00 



exists, yielding (|3.5|). We use the abbreviation 



e 2 r] 2 



2ir 2 



V 



-T 



in 



Let A denote the Lebesgue measure on 



I s and put 



Si :-- 



e 2 r) 2 
2^2" 



E 

(7=1,2 



T 



UJ 2 (U 2 + V 2 ) 



Xc dA. 



(5.6) 



Proposition 5.3. limjv->oo Si n = Si for each I G N. 
Proof. We will use dominated convergence. Put 



9i,N 

By the definition of V, we have 

Sin 



E 



e 2 7] 2 



2tt 2 



111 



^2 / 9i,Nd\. 

1 o ** 



(5.7) 



(5.8 



(7=1,2 



The functions y a , which contain the transverse vector fields, can be chosen to be continuous 
outside a closed Lebesgue null subset of M 3 containing the origin. Thus applying Lemma 



4.1 



(ii) on every factor on the right hand side of (5.7) one deduces that gi jy converges to 



the integrand in (5.6) almost everywhere. 
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Now each component of g^jy is bounded by 



E (vv^) 2 -^iAi< E (vv^) 4 ^iaj, 

its V A? itfc V ; W * 

where the inequality is obtained by applying Jensen's inequality to the convex function 



<fi : (0,oo) — >• M, x 1— > 1/y/x, yielding l/^/uJ] $S ( l/^/u)j for each j. Furthermore, we claim 
that 

(l7^) 1a^4-^1 Aj (5.9) 



Indeed, in the case G Aj, we replace Aj with the portion of the ball of radius V3/N that 
lies in the same sector. Then Aj is contained in that set, and 

,/Z (TJ75) =h~ ^ * = < if r < VS/tf. 

^ > 3 V 8 J 5 5-v/r 

whence the claim in this case. If Aj, we let (k/N,l/N,m/N) be the corner of Aj 
closest to the origin, and find 

/— =n < sup A .(l/VF) _ 1 y/(k + l) 2 /iV 2 + (Z + l) 2 /jV 2 + (m + l) 2 /^ 2 
V H /VW Jj ^ v^inf Aj (l/Vr) >/r ^/c 2 /iV 2 + Z 2 /7V 2 + m 2 /jV 2 

The last expression is maximal for (k, I, m) = (0, 0, 1) and takes the value v6 there. Thus 
{ 5T9J) holds. 

Hence every component of gi t N is bounded by (16/z^) 2 w -2 , which is locally integrable. 
This finishes the proof. □ 

A similar, but easier proof shows that 

lim £ V^f = p?- [ ~^—dX. (5.10) 



The integrals appearing in (5.10) and (5.6) can be calculated. We have 



I * C , dA = ^ln(l + f ). 
J uj 2 {lo + 77) c 3 



Together with the prefactor, and comparing with the definitions of sin tp, p and 7, we obtain 
the exponent in (3.5). For (5.6), note first that Yla=i 2 u & u a = ^3~jM^kk T follows from the 

_ T 

orthonormality of u%, 112 and k/\k\. Thus for each component of X)o-=i 2 I ZP^P^J 1 ) ^ 
the angular part can be calculated, and is found to be equal to Then 

E / W 2 ( "f; 2 } Xe dA = arctan (f ) J 3 . 

Taking the prefactor into account and computing the now trivial determinant gives the 
Z-th factor of the infinite product in (3.5). 

The final step is to justify the exchange of the infinite product with the limit A?" — > 00. 
To this end, we write, for each I £ N, 

2 2 
M t>N := \l z + Si jN , Mi := lim M^ N = \l 3 + S t . 
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Lemma 5.4. We have 

oo oo 

lim ff det(/ 3 + M hN ) = T\ det(/ 3 + Mj). 

2=1 2=1 

Proof. By the continuity of the exponential function, it will suffice to prove the result for 
the logarithms. Then the quantity of interest is given by 



In ( JJdet(I 3 + M l>N ) ) = £ln(det(I 3 + M hN )). 
\l=i / 

Since by continuity 



i=i 



lim ln(det(/ 3 + Mi N )) = ln(det(I 3 + Mi)) 

N— >oo 

for each / £ N, the result will follow by dominated convergence. 

At the end of the proof of Proposition 5.3 we have seen that Si t N is bounded by 
constant// 2 uniformly in N. Thus so is jv|| and hence its three eigenvalues. It follows 
that | det(J 3 + M^n)) — 1| ^ constant// 2 and therefore 

| ln(det(J 3 + M l>N ))\ < constant// 2 

uniformly in N . □ 



We have thus proved Theorem 3.2 



5.4. Removing the ultraviolet limit. We now investigate the limit 7 — > 00 in (3.5). 
Clearly, what we need to show is the convergence of 

00 

/( 7 ) := -2sinv9 pln(l + 7) + ^ln(l + fc,( 7 )), 

1=1 

where 

, / \ P 2 4 . flP\ 
'v7J = ~jj + ~ sm ^ arctan \~y) ■ 

We first note that 

0<fy( 7 ) < ^(l + ^7 sin ip) (5.11) 

/9 2 p 

< /i;( 7 ) /* hi := J2 + 2y sine/? monotonically for 7 00. (5-12) 



By (5.11), (^(7)) is summable for all 7, and thus 
/ 00 \ 00 

f{l) = E ^(7) - 2si n9 9 pln(l + 7) + ^ (in (l + fc,( 7 )) - fc,( 7 )) . (5.13) 



\2=1 / 2=1 

We now use the formula 

00 

C = lim 



^aretanft-ln. 

7T ^ / / 

2=1 



(5.14) 



for the Euler/Mascheroni constant C. The first bracket above is then equal to 

00 2 / 2 00 1 \ 

E ^ + 2s[ni P P lnp + ln + y arctan(^) - ln(p7) , 

2=1 V 71 2=1 / 
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and by ( |5.14[ ) converges to ^- + 2sin(</?)/9(C+ln(p)). The last sum in (5.13) converges 

to Y2i=i ( m (l + ^i) ~~ ^i) ^ (5.12) and monotone convergence, since ln(l + x 
monotone decreasing. Combining these two results, we obtain 



x is 



/:= lim f(j) = 2simp p(C + ln(p)) + V(ln(l + /iA -2^ simp 

■y— >oo ^—^ V v 'I 

= 2simp p(C + ln(p)) + In I J|(l + /i;)e~ 2 i simp J . 



Taking into account that 1 + hi = (1 + %ie lLp )(1 + jie lLp ), and writing 2isin</? = e 
e~ lip , we find 



^2p ln(p) sin </3 



pe 



Cp i e 



' n i 1 + 



-i e 



c i ' 



Using the representation 



Hz)- 1 = ze Cz 



z=i 



1 e ; 



for Re z > 0, 



for the Gamma function on the right half-plane, we arrive at (3.6). 

6. The effective measure of states 
In this section, we investigate the inverse Laplace transform of Z((3). As a first step, 



we prove the second part of Theorem 3.3 



6.1. Complete monotonicity. By Bernstein's theorem, a function / on [0, 00 [ is the 

Laplace transform of a positive Borel measure if and only if it is completely monotone 
(cm.), i.e. (— l) n 9™/(x) ^ for all x 0. To show that Z(/3) is cm., we use (3.6) in 
order to get 

1 



InZ(P) = Inp- ln(2vr) - 2 sin ip plnp + lnT(pi e~ lip ) + \nT(pi e"^ ). 



We now use Binet's formula 
lnT(z) = 



00 — tz 



t 



1 1\ ln(2vr) . 1 \ i / \ 

---) dt+^^ + (z-l)ln(z)-z, 



valid for Re (z) > 0. Since ln(z) = ba.(z) for all z £ C, we obtain 



with 



9(t) :-- 



Re 



lnZ(/3) 

1 1 

~ r ~ 2 



e~ tp dt - %p p 



(6.1) 



7T 

e~ tlf t, ttp := 6(sm</? + (— — ip) cosip). 



t := 1 



Below we will show that g(i) for all t 0. This immediately implies that the first 
term on the right hand side of ( |6.1[ ) is cm. Since the product of two cm. functions is cm. 
also exp(/) is cm. if / is. Thus Z(f3) is cm. as the product of exp(/ °° e~ tp g(t) dt) and 
the the clearly cm. function exp(— t v p). The claim g(t) ^ follows from 
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Lemma 6.1. For all w G C with Re (w) > we have 

1 1 



h(w) := Re 
Equality holds if an only if Re (w) = 0. 



1 



w 



> 0. 



Proof. We write w = u + iv. If u = 0, then h(iv) 
calculation yields 



1 — COS V 

1—2 cos v+1 



0. For u > 0, a direct 



(u 2 + u 2 )(l 



-2u 



2u(l + e" 



-2u 



2e n cost;) 



2(1 - 2e"" cost; + e" 2 " )(u 2 + u 2 ) 
(u 2 + v 2 ) sinh(u) — 2«(cosh(it) — cos(u)) 
2(u 2 + f 2 )(cosh(-u) — cos(v)) 



with nonnegative denominator. We need to show that the numerator r(u, v) is positive if 
u > 0. By symmetry it suffices to treat v 0. Take first t> = 0. Then d u (r(u,0)/u) = 
wcosh u— sinhu, which at u = equals 0. Since <9 2 (r(u, 0)/u) = u sinhu > 0, u \— > r(u, 0) > 
for u > 0. Thus the result holds if d v r(u, v) = 2(t>sinh-u — lisin-u) is nonnegative for all 
u > 0, v > 0. But this is clearly true, since (sinh u)/ii > 1 > (sinv)/v for all u, v ^ 0. □ 



6.2. Analytic continuation. Here we prove Theorem 3.5 Let us define 

Y(p) := Z(P)e^. 

Since p = jfe, the inverse Laplace transforms of Y(p) and Z(/3) are related by translation 
and scaling. The key observation is that, due to (6.1), Y is itself the exponential of a 
Laplace transform, 



oo 1 

Y(p) = exp£g(p) = 1 + £ -£<?*». 

n! 



n=l 



Above, £ denotes the Laplace transformation, and 



h* k(x) 



h(t)k(x - t) dt 



(6.2) 



(6.3) 



is the convolution of functions supported on [0, oof. Since g is bounded on compact in- 
tervals and nonnegative, the final sum in ( |6.2| ) converges uniformly on compact intervals. 
Let us define 

oo 1 



71=1 



By Lemma 6.1, all terms of the above sum are positive. Thus £g = Y — 1 holds by 
monotone convergence, and the uniqueness of Laplace transforms shows £ _1 1" = 5o + g. 
Moreover, a direct calculation shows that g(t) — > ^ sim p as t — > 0. Thus the same holds 
for g, and we have shown the first part of Theorem 3.5, with <p{oj) = ^p(^j-(uJ — ^v)) ^ or 
oj > uj^p and zero else. 

We now investigate the analytic continuation of g. First note that the analytic contin- 
uation of g is given by 



9(z) 



1 



+ 



2 sin <p 



1 — exp(— ie li P z) 1 — exp(i e llfi z) 



1 
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Lemma 6.2. g is a meromorphic function on C with simple poles at the zeros of the 
denominators 

q j = 2-Kje~ i ' p and q] = 2irj e iip , j € Z \ {0}. 
The residues are Res(g; g, ) = Res(g;qj) = =|. 

The proof is routine. Note only that we have seen already above that the common zero 
z = of the denominators actually is a regular point of g. 

Let us now consider the analytic continuation of the convolutions. Define V := {z : 
z = qj or z = qj for all j £ Z \ {0}} and 5 := {z : z = se^ for real s with |s| > 2ir}. 
For brevity use A(S,V) to denote the set of functions / that are analytic on C \ S such 
that the components of S \ V are cuts for / emanating from the points of V . Analytic 
functions on C \ V are elements of A(S, V). — For any function / and any subset K of C 
let || /|| k denote the supremum of |/| on K. 

Lemma 6.3. Let h £ A(S,V) and let k be analytic on C \ V . For z£C\5 set 

F(z) := f h(0K* -C)dC 



o 



integrating along the straight line joining and z. Then F £ A{S,T') and F is an analytic 



continuation of h * k as defined in (6.3). Let K be any compact subset ofC\S and let 
K denote the union of all straight lines joining the origin with some point of K. Then 
\\F\\k < \\id\\ K \\h\\x \\k\\ft. 

Proof. Plainly, F agrees with h * k on the real axis and is differentiable at all z E C \ S. 
Hence F is an analytic continuation of h*k on C\S. — We turn to the analytic continuation 
F of F, e.g., from above across the cut between the points qj and qj + i with j £ N. Let 
z £ Si := {z : z = se"^ for s > 2ir} be between qj and qj+i- Join the points and z 
along S\ but avoiding qi, respectively z — qi, for / = 1, . . . , j, by making a small detour 
above, respectively below, S\. By analogous curves ^ z one joins to z for every z in the 
open disk D centered at \{qj + qj+i) with radius it. Let h be an analytic continuation of 
h from above Si across the cuts along Si. Then 



F(z) := / h(()k(z -()d( 

A* 

defines an analytic function on D. It extends F, since for z £ D above Si the closed 
curve composed by 7 Z and the straight line from z to does not contain any singularity 
of C — > h(C)k(z — C)- — The remainder is obvious. □ 

By the lemma, g* n £ A(S,V) for all n £ N. Moreover, \\g* n \\ K < NHk" 1 llslll- 
Obviously, a similar estimate holds more gene rally for the analytic continuations F of F 



on compact K C domF. Hence the series in (6.4) converges uniformly on compact sets 



implying that the limiting function belongs to A(S,V). This proves the second part of 



Theorem 3.5 when taking into account that the translation and scaling takes qj into pj. 



Let us comment on the cuts of F from Lemma 6.3 Even if h and k are analytic on 
C\V with poles at points of V, the convolution F may have cuts. More precisely, e.g., if 
z lies on the cut between qj and qj+i with j £ N then 

3 

F+(z) - F-(z) = -2ir[Y,[k(z ~ qi) Res(h; qi ) + h(z - qi) Res(k; qi )} 

i=i 
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where F + {z) and F_{z) denote the limit values of F at z approaching z from above and 
from below the cut, respectively. This is an immediate consequence of Residue Theorem 
integrating the meromorphic function £ — > M z (() := h{C,)k{z — £) along the simply closed 
curve, which is symmetric with respect to S\ and which joins to z by j z . — In case of 
F = g * g the above formula yields the jump function 



1=1 



6.3. Analysis of the singularities. Here we prove Theorem 3.6 Again it suffices to 



analyse g (6.4) instead of <f> as <j){uS) = ^^(^(w — oj v )) for uj > uo^. We refer to the 



analytic continuations of the convolutions g* n and of g defined in Section |6.2 
First note the following formula for z = se lx with s > and |x| < ip 



1 



1 



o C - Qj z ~ C - Qk 



dC 



-27ri + ln(z - qj) - lngj + ln(z - q k ) - \aq k 
Qj ~qk z - qj - q k 



1 + 1 



which follows from the partial fraction expansion -A i — = — , , 

evaluating the primitives ln(£ — qj) and — ln(z — £ — q k ) of (C — qj)~ l and (z — £ — qk)' 1 
respectively. Note that the second term in (6.5) is regular at qj + q k = qj+k- Next let us 
define recursively the coefficients 



j'-i 



Cjl — Cjl — ., Cj,n+1 — ^ ^ CklCj—k,n> Cj >n ^-% — ^ 1 CklCj—k,n (6-6) 

3 k=l k=l 

for all j, n E N, where the void sums for j = 1 are zero. Set 

Cjn ~ / \ Ci 



Sjn(z) :- 



2iri(z — qj) ' 



l -jit 



2iri(z — qj) 



In the following, we will concentrate on the forth quadrant of C. Subsequently it will 
be easy to extend the result to the right half-plane. Fix N £ N and let z = se lx with 
< s < 2ttN and -f < % < 0. 



Proposition 6.4. Then 



N 



g* n (z) = ^ s jn{z) + L n (z) for \x\ < V, 
3=1 

N 

g* n (z) = ^ Sjn(z) + L n (z) for \x\ > <p 
hold with \L n (z)\ and \L n (z)\ bounded by A n (l + I In sin \ ip — \x\\\ S ) f or some constant A. 



Proof. We proceed by induction. The statement for n = 1 follows from Lemma 6.2; indeed 



subtracting from g the first order poles at qx, . . . , qw leaves a bounded function L\. Let us 
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9 



*(n+l) 



now assume that g* n has the asserted decomposition and consider the case \x\ < f- Then 

( z ) = Y] / Sji(C)s kn (z - C)d( 
3,k=i Jo 

+ E / (MOLniz - C) + SjniQLtiz - 0) dC 



(6.7) 



By (6.5) we find 



Sji(()s kn (z - C) dC 



2vri(z 



+ Ljk n (z), 



(6.8) 



where i^-fc n is regular at qj +k - The first term above contributes to Sj +k ,n+i- We will show 
below that none of the remaining terms entering g*( ra+1 ) has any first order poles, and 
thus by collecting all terms with j + k = m we obtain the recursive equation for c^n+l- 
The calculation for §j n , i.e. for \x\ > f, is very similar. The only difference is that the 
residue of the pole at qj+k in (6.5) is 2tt\ instead of — 27ri, a difference which is due to two 
jumps of 2tt\ of the logarithmic terms at the cut along the negative real axis. This gives 
the additional minus sign in the recursion for Cj n . 

We turn to L^ n . As mentioned above, there is no singularity at qj+k- There are two 
logarithmic singularities at qj and g&. Other than that, Lj^ n is bounded. Set 8 := \ip— \x\\ 
and h(5) := |ln(sin<5)|. Then \Ljk n (z)\ ^ Cj\Ck n K(l + h(5)) for some constant K. It 
is immediate from (6.6) that Cj n = for j < n. Thus there exists some constant B, 
independent of n, such that 

N 



£ \Ljkn{ z )\ ^B(l + h(S)). 



j,k=l 



Now we tackle the second line of (6.7). By the induction hypothesis, |L n (z)| ^ A n (l + 
h(5) s ). Thus 



Ijn(z) :-- 



Sjl (0L n (z-()d( 



< A n h(sy 



Sjl (ye^)\h(8)^dy + A n / \s n (ye^)\dy. 
o Jo 



For s < 7T, the integrands on the right hand side above are bounded, and hence Ij n (z) K 
DA n (l + h(5) s ) for some constant D. For s > it, we decompose the domain of integration 
into y < it and ir < y < s. On the first interval, the integrands are bounded with the same 
result as above. On the second interval, we replace h(5)~ y by h(8)~ v . The integral over 
s\j alone is clearly bounded by ^(l + h(5)) for some constant E, and we estimate crudely 

N 

Ijn(z) < 3NA n (D + E)(l + h(5) s ). 

3=1 

Finally, since L\ is bounded, we have 

N 



E 

3=1 



s in (C)Li(z-C)dC 
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where the constant F does not depend on n. Altogether, we find |iv n +i(z)| 3NA n (D + 
E + B + F)(l + h(5) s ). Clearly, setting A = 3N(1 + D + E + B + F), this is bounded by 
A n+1 (l + h(5) s ). □ 

In order to extend this result to the right half-plane one has to take account of the 
poles g„- of g, too. This amounts in replacing Sj n (z) by Sj n (z) + Sj n (z). Let M n denote the 
remainder in place of L n . It satisfies the same kind of estimate with some new constant 
A. Using (6.4) and the fact that Cj n = for j < n, we now have 



N 



j=l \n=l 



ill 



1 



+ 



1 



+ M(*) 



27ri(,z — 2ir\(z — qj) 

with M := ^~=i and \ M ( Z )\ < e A (1 + /i(5) s ) for -99 < x < <£■ The same formula 

with dj n and M holds for \x\ > ip. Obviously h(5) s can be replaced by | ln5\ 2nN . 

Proposition 6.5. Define Cj n and Cj n as in ( |6.6[ ). Then, for all j 6 N, 

3 / • , o\ j 



n=l 



i + 2 
2 



and 



n=l 



J" 



-1 



i/j ^ 3, and otherwise. 



Proof. We introduce the generating functions -F n (x) = Sfci^'n^" 7 - Then, iq( 
Sj^=i j 2 ^ = — 31n(l — x), and F n = i*Y\ The last statement follows from 



Fi(x)F n (x) 



00 /i-l 
i=i Vfe=i 



3=1 



F n+1 (x). 



Now recall that c, n = for j < n. Then 



1 



n=l n=l 



n! 



iai ( e F ^ - 1 



x=0 Jl 



x=0 



l^ e -31n(l-x.) 



j\ x {l-xf 



x=0 



i + 2 



For the generating function F n of (cj n ), a similar calculation leads to F n = (— l) n+1 .F ] n 
with F\ = F\. Thus, as above, 



00 „ 

ECjn_ 



1 



gj e 3\n(l-x) 



n=l 



x=0 



r^(l-x) ; 



x=0 



This equals 3 for j = 1, —3 for j = 2, 1 for j = 3, and zero otherwise, as was claimed. □ 



The stated analyticity properties of h^{z) and h^{z) from Theorem 3.6 are true by 
the fact that both <f> (cf. Theorem 3.5) and £j possess them. This concludes the proof of 
Theorem 13.61 



7. Appendix 

Here we show that Definition [2] of a resonance is equivalent to that given e.g. in |X6|, 
XII. 6]. We start with a Lemma that connects the analyticity of a measure's density with 
properties of its Stieltjes transform. 
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Let \x be a finite Borel measure on [0, oo[ and / := /x its Stieltjes transform. Note that 
/ is holomorphic on C \ [0,oo[ and satisfies f(z) = f(z). We recall the classical inversion 
formula valid for t > 0: 

1 f t+s 

/x([0, i]) = lim 5 ^ 0+ lim e ^ + / Im f(s + ie) ds. (7.1) 

Lemma 7.1. Let to > and C7 C C an open disc centered at to with radius smaller than 
to- Then the following two statements are equivalent: 

(i) there is an holomorphic function F on U which equals f on {z £ U : Imz > 0}. 

(ii) fi on U n M is absolutely continuous with respect to Lebesgue measure, and its 
density is the restriction on U n R o/ a holomorphic function (j) on U . 

If (i) and (ii) hold, then obviously f(t) : = lim e _ > o+ converges uniformly on compact 

subsets ofUnM, and F = f + 2iri(J) holds on {z £ U : Im 2 < 0}. 



Proof. Assume first that F exists. Then (7.1 ) yields immediately for t%, ti in J7nR, t\ < t2~. 



M]*i)*2]) = \ // 2 ImF(s)ds. This implies that /i on U D R is absolutely continuous with 
respect to Lebesgue measure and that the density is given by t -)■ ilmF(t). Then 
0(z) := {F{z) — F(z)) is its analytic continuation on J7. For z £ U with Im z < one 
has F(z) = /(i) = f(z), whence = /(*) + 2vri0(z). 

Now assume the existence of <f>. For z £ U set := /(z) if Imz > and -F(z) := 

f(z) + 2yri0(z) if Imz < 0. Then F is holomorphic on U \ R. Fix i £ [/ fl 1. Let 5 > 
such that {t — 5, t + 5} C C7. We define three paths in C. First 71 (s) := s for s E [0, 00 [. 
Then 72 differs from 71 only in that it joins the point t — 5 to t + 5 not by the straight 
line but by the semi-circle through t + i5. Finally the closed path 73 joins t — 5 to t + 5 
forwards by the straight line and backwards by the semi-circle through t + i5. Then for 
< e < 5 one has 

dfj,(z) f d[x(z) f (p(z)dz 



, 7l z - (t + ie) J l2 z - (t + ie) 
by the residue theorem. Therefore 



L 3 z-{t + ie) 



lim £ _> + f{t + ie) = / + 2vri 



exists. Similarly lim e _j.o+ f{t ~ i £ ) is shown to exist. 

Set F(t) := lim e _ > o+/(i + ie)- Thus F is defined on the whole of U. It remains 
to show that F stays holomorphic. By the following lemma F(t + ie) — F(t — ie) = 
f(t + ie) — f(t — ie) — 27ri </>(t) + 27ri (</>(i) — <^>(i — ie)) — >■ as e — >■ 0+ uniformly on compact 
subsets of U H M. From this the premises on F of Morera's theorem easily follow, whence 
the result. □ 

Now the equivale nce o f Definition [2] with the traditional definition of resonances follows 

Mo 



by applying Lemma 7.1 to [i = and \i = //q : choosing U such that the density 4> of 



or is analytic on U, we find that the continuation F on U of the resolvent to the 
second Riemann sheet is given by / + 2-7ri0. 



We close with a lemma showing that a strong version of (7.1) holds if the density of \i 
is continuously differentiable. 
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Lemma 7.2. Let fx be absolutely continuous with respect to Lebesgue measure on ]A,B[ 
for < A < B and let the density <p be continuously differentiate. Then <p(t) = lim e ^o+ 
^lmf(t + ie) holds uniformly on every compact subset of]A,B[. 

Proof Put x(s;t,e) := • Then £lm/(t + ie) = Jj 0oo[ x(s; t, e)d/i(s). Let A < 

Ai < Bi < B and t € [A ± ,Bi]. 

i) First j [B oo[ x(s;t,e)dfi(s) < ^-BxY Jjo,oo[ d M s ) uniformly as e ->• 0+. Similarly 
this holds for Jj Q ^ %(s; t, e)d/x(s). 

ii) Next 

1> / XI s ; e ) ds = — (arctan arctan )> 

Va tt e e 

\, B-Bx At -A, 

^ —(arctan arctan ) — > 1 

7T e e 

as e — )■ 0+. This implies that ^ e ) ds tends uniformly to 1 as e 0+. 

iii) Because of i) it remains to show that \ c/)(s)x(s;t, e)ds — 4>(t)\ < \J^(4>(s) — 

(j>(t))x(s; t, e) ds\ + (f)(t)\ x( s ; t, e)ds — 1| tends uniformly to as e — >■ 0+. This holds 
true for the second summand because of ii). As to the first summand the mean value 
theorem yields \4>(s) — 4>(t)\ < c\s — t\ with c := sup{0'(r) : r € [^4i,5i]}. This finishes 
the proof since \s - t\x(s;t,e) ds = | f^~* \r\(r 2 + e 2 ) _1 dr < ^ f^r(r 2 + e 2 ) _1 dr = 
L In -»• as e -> 0+. □ 
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